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We present the first realization of Stochastic TDHF, a the- 
ory which goes beyond pure mean-field dynamics, embracing 
dissipation as well as fluctuations. Applications to heavy-ion 
collisions in the Fermi energy domain are given and analyzed 
in terms of dissipative features and Intermediate Mass Frag- 
ment production. 



Since its early introduction in quantum dynamics 
Time-Dependent Hartree-Fock (TDHF) is the basic 
selfconsistent theory for many dynamical processes in 
physics and chemistry. In nuclear physics, it has found 
since long widespread and fruitful applications in collec- 
tive motion and heavy- ion reactions 0-0]. Recently, a 
new field of applications to the nonlinear electron dynam- 
ics in metal clusters has opened up ^-|| which brings up 
challenging questions in connection with new experimen- 
tal investigations 

Mere TDHF is nevertheless too restricted in its 
degrees-of-freedom if one goes beyond the regime of low- 
energy (collective) dynamics. For example, TDHF be- 
comes insufficient for the description of central heavy-ion 
collisions in the Fermi energy domain But it is just 
these collisions which have unravelled a bunch of exciting 
phenomena, particularly from the generation of experi- 
ments using 47T multidetectors |jq ], and there exists up 
to now no theoretical tool which is definitively validated 
in this field (16 1 . 

Extensions of TDHF to include dynamical corelations 
are formally and computationally feasible |l4|Jl7| , p^ | . But 
these extensions still remain restricted to the low-energy 
regime, where only few channels are active in the dy- 
namics. Statistical approaches are more appropriate for 
higher energies and there exist attempts to complement 
the TDHF equation by a collision term |l^,^(|. This 
raises, however, considerable technical difficulties due to 
the large number of rarely occupied orbitals needed in 
such an approach. Moreover, such a collisionally ex- 
tended TDHF is still bound to one mean path and grossly 
underestimates the large statistical fluctuations of the 
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mean field. 

The semiclassical variant of collisional TDHF, the 
Boltzmann - Uehling - Uhlenbeck (BUU) approach, has 
been more successful in describing average trends in 
medium-energy heavy ion collisions fl2l]]22|| . But it suf- 
fers also from insufficiencies, particularly when consider- 
ing experiments in which fluctuations are playing a cru- 
cial role such as in the formation of Intermediate Mass 
Fragments (IMF's) jl5| or in the production of rare parti- 
cles well below threshold p3fl . The Boltzmann-Langevin 
equation (BLE) has thus been introduced as an extension 
of BUU which incorporates (possibly large) dynamical 
fluctuations by switching to a full ensemble description. 
It should allow a pertinent description of variances and 
widths of observables 24 1. But huge computational (and 
even some formal) difficulties hinder up to now large scale 
applications of these BLE approaches. 

A related theory for generating and propagating an en- 
semble of quantum states is Stochastic Time-Dependent 
Hartree-Fock (STDHF). It goes back to the basic quantal 
mean field dynamics and treats dynamical correlations in 
Markovian approximation by occasional jumps from one 
TDHF trajectory to another one, thus generating an en- 
semble of TDHF trajectories in the course of time p5| . 
It is interesting to note that similar trajectory hopping 
techniques have also succesfully been applied to electron 
dynamics in small molecules 1 26 . From a formal point of 
view, STDHF is the basic theory from which collisional 
TDHF, BUU and BLE can be derived by appropriate 
subsamplings p5[ . And technically, it offers a particu- 
larly efficient scheme to evaluate at once collision as well 
as fluctuation terms p4| by means of trajectory hopping. 
The ensemble of pure Slater states is a more efficient 
way of handling thermalization because it circumvents 
the complication to carry an extra load of nearly unused 
single particle wavefunctions. Furthermore, its quantal 
nature avoids the many (formal and numerical) difficul- 
ties encountered in a proper justification [p7|,[16[ and sim- 
ulation of semi-classical kinetic equations |p8| in the con- 
text of heavy-ion collisions. 

The key issue in the STDHF scheme is an efficient 
selection of the jumps. Strategies have been proposed, 
based on the principle of maximum overlap with the 
corresponding correlated state and evaluated practically 
with linear response techniques [^,|30|]. But before at- 
tacking the problem with these elaborate (and costly) 
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optimization principles, it is worthwhile to investigate the 
capabilities of STDHF as such, by using an intuitive ap- 
proach to the jumps. Our procedure is inspired by BUU 
and BLE, where the (hardly tractable) 5-fold local Boltz- 
mann collision integral is replaced in a robust manner by 
sampled collisions of test particles, simply parametrized 
by one parameter, a total effective nucleon nucleon cross 
section a [plf . It is the aim of the present letter to pro- 
pose such a simple, intuitive, and efficient scheme for 
generating the ensemble of STDHF, and to apply it to a 
realistic test case, collisions of 16 on 16 at medium en- 
ergies, where pure TDHF only shows transparency, while 
dissipation should lead to incomplete fusion, fission or 
fragmentation ]l3[ |. 

First, let us recall briefly the evaluation of the collision 
term in BUU (2^j3l[] . The total cross section a is rescaled 
to an effective cross section for the test particles <7 e ff = 
va, where v is the number of particles per testparticle 
PjJ . For all pairs of test particles "(1,2)" one checks 
whether " 1" and " 2" will ap proach each other closer than 
a scattering length d a = ^/ceff/i") within the next time 
step. If yes, "1" and "2" collide, i.e. their momenta 
are changed (fix, Pi) — > (p'iiP^i while preserving total 
momentum and total kinetic energy. Finally, the Pauli- 
blocking factors (1 — /{)(1 — f-i) are computed at the new 
phase-space points (r\,Pi) and (r2,p 2 ) an d sampled by 
Monte-Carlo techniques. 

A TDHF state |$) consists of single particle wavefunc- 
tions ip a (r, t) where each one represents one physical par- 
ticle and is spread over all space, carrying a local density 
distribution p a (r, t) and a local momentum distribution 
p a (r,t) = j a (r,t)/p a (r,t). Collisions, on the other hand, 
are local. We thus split our space grid into "test volumes" 
V W which are going to play the role of the test particles 
in BUU. A spatial filtering function (r) is associated 
with this test volume and we define V^-test volume in- 
tegration as J v(i) dr... = J drC/W(r).... For a given VW, 

we hence evaluate the "local" average positions r a — 
J v(0 rp a dr/N a l > and momenta p$ = / y(i) j a (v)dv/N a l) , 

where N a = J*y(i) p a dv. These positions and momenta 
are then handled in exactly the same manner as in BUU. 
The only modification is that the particle content v now 
becomes a dynamical quantity to be evaluated as the av- 
erage fraction of the two scattering particles a and (3 in 
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We thus obtain a jump of the associated momenta 



local momentum and kinitic energy as 



obeying local conservation of 
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We need now to modify the phase profile of the single 
particle wavefunctions <p a and (pp such that local average 

over produces the new momenta (p'a\pg^)- This is 
achieved by the transformation 



tf/ a (r,t) = exp(-itf a -p a )rg®(rj)tp a {r,t) (3) 

and similarly for ipp. 

As in BUU, the collision may end up in a region of 
phase space where particles are already present. This 
is measured by the overlap of the old $} and new |$') 
TDHF states. The Pauli blocking factor is the amount 
of still available space, which is then simply 1 — |($'|$)| 2 . 
This blocking factor is sampled by Monte-Carlo tech- 
niques. In contrary to the classical BUU scheme we need 
here two more steps for a final clean-up. First, the new 
single particle wavefunctions (p' a and tp'^ are each orthog- 
onalizcd on all other ipy. And second, we need to restore 
the original value of the total energy because we find a 
slight modification of the total energy (although most of 
the energy matching has been accounted for by conserv- 
ing the " classical" kinetic energy in the collision) . This is 
achieved by energy rescaling with an iterative imaginary 
time step |$') — > (1 — Sh) n \&') where h is the actual 
mean field Hamiltonian and 5 a (small) shift parameter 
to be chosen such that total energy is recovered in a few 
iterations n ~ 4 — 5. This altogether accomplishes the 
stochastic jump of the STDHF trajectories. 

The scheme can be simplified for the initial stages of 
a collision of two small nuclei. The single particle wave- 
functions in each one of the colliding nuclei are initially 
well localized, such that we can skip the subdivision into 
test volumes and treat the collision with respect to the 
straightforward averages f a and p a of the single parti- 
cle wavefunctions. This amounts to letting the space- 
filtering factor C/W — > 1 and the trivial renormaliza- 
tion factor v = (Na ^N^) 1 / 2 — > 1. This simplification 
is surely applicable for the initial stages of the collision 
where most of the statistical ensemble is built up and 
where most dissipation takes place It may overesti- 
mate collision rates after the compound stage. But this 
is mostly a quantitative problem. The description will 
remain pertinent at least at a qualitative level. 

For a first exploration, we have implemented the 
above scheme in its simplified version (using global av- 
erages rather than local ones) into a 3D nuclear TDHF 
code without any symmetry restriction. The mean field 
Hamiltonian is generated from a Skyrme force with 
a density-dependent zero-range interaction, as used in 
BUU studies pq| . As a test case, we consider 16 O on 16 O 
collisions. 

Figure |l| shows the results for a head-on collision 
16 0+ 16 at 5QMeV/u beam energy. The STDHF 
scheme has been carried up to an ensemble of 100 tra- 
jectories. The upper panel shows the time evolution of 
the ensemble average momentum quadrupole moment 
Q20 = (k 2 Y2o)- It is known to provide a robust indi- 
cator of the dissipative phase leading to the formation of 
an excited compound nucleus (distinguished by a small 
deformation) j32U23] . A TDHF and a BUU result are also 



plotted for comparison. The TDHF calculation leads in 
this case to a perfect transparency which is reflected by 
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typical mean-field oscillations of Q2o- In contrast, the 
STDHF calculation, as BUU, leads to the formation of 
a compound nucleus as can be seen from the damping 
of Q20 to 0. The damping rate t c in STDHF, how- 
ever, is somewhat larger than in BUU. This rate may 
be further analysed by considering the number of succes- 
ful two-body collisions function of time, see lower 

panel of Figure [jj Two-body collisions are peaked around 
t = 30fm/c, which roughly corresponds to the time of 
maximum overlap of the two colliding nuclei. The BUU 
calculation provides a smaller peak value of n c , but a 
broader peak, such that, after all, the integrated number 
of two-body collisions during the overlap phase is com- 
parable to the result of STDHF. 
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It is interesting to note here the difference between 
BUU and STDHF, before overlap. While in STDHF Pauli 
blocking is perfectly effective, leading to exactly no two 
body collisions (n c = 0), there remains a residual back- 
ground of 2-body collisions in the BUU case, due to the 
approximate treatment of Pauli blocking and to the fact 
that the test particle sampling only provides approximate 
ground states, even with large numbers of test particles 
per nucleon. Finally, the middle panel of Figure [l] dis- 
plays the variance 020 associated to the quadrupolc mo- 
ment Q20J as calculated from the ensemble of computed 
events. Note that a 20 exhibits a marked bump at the 
time of maximum overlap. Such a bump in variance is 
characteristic for a transient regime as demonstrated in 
BLE-based calculations |J|||. 

As stressed above, a major advantage of STDHF is the 
fact that it provides at once an ensemble of TDHF tra- 
jectories, which contain the proper (possibly large) fluc- 
tuations of the mean-field. As a first illustration of the 
capabilities of STDHF, we analyze the emerging ensem- 
ble in terms of the formation of IMF's. The IMF identifi- 
cation is done with a percolation-type algorithm applied 
to the density of matter inside the computing box, as in 
p4j . This analysis allows to evaluate, mass, charge, ve- 
locity,... of the IMF's. The upper panel of Figure |^ shows 
the fraction /1 of events leading to incomplete fusion (1 
fragment) as function of time. Results from four different 
16 O + 16 O collisions at beam energies 50 and 75 MeV/u 
and impact parameters 6 = and 2 fm are shown. Just af- 
ter touching (t~2fm/c), fi goes to 1 which corresponds 
to the overlap phase. In a later stage (t>40 — 50fm/c) 
the system evolves to either incomplete fusion (/1 = 1), 
fission (2 IMF's), or fragmentation (3 or more IMF's). 
At 50 MeV/u, both impact parameters lead to almost 
100% incomplete fusion (average mass < A >~ 22 ±1 
and charge < Z >~ 11 ±1). Note that isospin symme- 
try (N=Z) is preserved here, without being enforced. At 
75 MeV/u, fission and even fragmentation mechanisms 
take over. The detailed division into various fragments is 
shown for the most " fragmenting" case in the lower panel 
of Figure ||. Fission seems still to dominate. But there 
is a sizeable amount of events with three fragments and 
occasionally more than that. 

It should finally be noted that the computational ef- 
fort involved in one STDHF event is comparable to the 
effort required for performing an acceptable BUU calcu- 
lation. This has to be put in perspective with the differ- 
ences in the physical contents of the 2 approaches (quan- 
tal versus semi-classical, full fluctuations versus average 
description) . 



FIG. 1. Time evolution of a collision lb O+ lb O with 
impact parameter b — Ofm at 50MeV/u beam energyi 
(STDHF: ensemble of 100 events, BUU and TDHF: 1 (aver- 
age) calculation). Uppermost block: the quadrupole moment 
Q20 = (k 2 Y2o)- Second block: variance of Q20 in the STDHF 
ensemble. Third block: number of successful 2-body collisions 
per fm/c in STDHF and BUU. 
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of his 3-D TDLDA code and P. Bonche is acknowledeged 
for discussions in the final stages of the work. 
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FIG. 2. Time evolution of the fraction of events 
(100 events ensembles for each reaction) leading to incom- 
plete fusion (upper panel) and to fission of fragmentation 
(lower panel) for 16 0+ ie O collisions with impact parameters 
6 = and 2 fm at 50 and 75 MeV / u beam energy. 

We have presented a first application of STDHF as a 
propagation scheme which is capable of describing colli- 
sions and fluctuations beyond a mean-field approach by 
generating an ensemble of TDHF trajectories through 
occasinal jumps. An intuitively motivated and well man- 
ageable jump algorithm was developed in close analogy 
to the semiclassical cousin of collisional TDHF, the BUU 
equation. The first results for collisions 16 0+ 16 are 
very encouraging. We see a physically reasonable damp- 
ing of collective motion leading to formation of a com- 
pound nucleus. This proves that STDHF is an efficient 
way of sampling a collision term in the case of quan- 
tal propagation. Beyond that, STDHF generates at the 
same time an ensemble of trajectories with the capabil- 
ity of embracing large mean-field fluctuations, which cor- 
responds to the Boltzmann-Langevin treatment in the 
semiclassical domain. This also allows a reliable descrip- 
tion also of variances. The results of our test case show 
nicely the growth of the variance related to strong dissi- 
pation. Moreover, these results have been obtained with 
much less expense but more physics than needed for a 
comparable BLE calculation. 



[1 
[2 
[3 

[4; 

[5 
[6 

ir 

[8 
[9 
[10 

[11 
[12 

[13 

[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22 
[23 
[24 
[25 
[26 
[27 
[28 

[29 
[30 

[31 
[32 
[33 
[34 



P. A. M. Dirac, Proc. Camb. Phil. Soc. 26 (1930) 376 

J.W. Negele, Rev. Mod. Phys. 54(1982)912 

K. Goeke, P.-G. Reinhard (Eds.), TDHF and beyond, 

Lect. Notes in Phys. 171, Springer, Berlin, 1982 

D. Klakow et al, Z.Phys. A351(1 995)391 

F. Bertsch et al, preprint lanl: frucl-th/9692050 



Calvayrac et al, Phys.Rev. B15(1995)rl7056 
F. Bertsch, ISSPIC8, july 1996, Zeit. Phys. D, in press 
Calvayrac et al, Ann. Phys. (NY) 255 (1997) 125 
Chandezon et al, Phys. Rev. Lett. 74, 3784 (1995) 
Reiners et al, Zeit. Phys. D, in press 
M. Barat et al, Zeit. Phys. D, in press 
K. Wohrer et al, J. Phys. B29(1996)L755 
K.T.R. Davies et al, in Treatise of Heavy-ion science, 
vol3, Ed. D. A. Bromley, Plenum, 1985 
P.-G. Reinhard et al, Nucl. Phys. A398(1983)141 

B. Tamain et al, Les Houches, 1996, LPCC 96-16, in press 
E. Suraud, Ann. Phys. (Fr.) 21(1996)461 

R. Balian, M. Veneroni, Ann. Phys. (NY) 135 (1981) 270 
J.M. Marston et al, Phys. Rev. Lett. 54 (1985) 1139 

C. Y: Wong, K.T.R. Davies, Phys.Rev. C28 (1983) 2321 
H.S. Kohler, B.S. Nilsson, NUcl. Phys. A477 (1988) 190 
G.F. Bertsch, S. Das Gupta, Phys. Rep. 160(1988)190 
W. Cassing et al, Prog. Part. Nucl. Phys. 25(1990)235 
J. Julien et al, Phys. Lett. B264(1991)269 

Y. Abe et al, Phys. Rep. 275(1996)49 

P.-G. Reinhard et al, Ann.Phys. (N.Y.) 216(1992)98 

J.C. Tully, J.Chem.Phys. 93 (1990) 1061 

P. L'Epplatenier et al, Ann. Phys. (NY) 244(1995)426 

P.-G. Reinhard, E. Suraud, Ann.Phys. (N.Y.) 239 (1995) 

193; and 239 (1995) 216 

P.-G. Reinhard, E. Suraud, Nucl. Phys. A545(1992)59c 
P.-G. Reinhard, E. Suraud, Heavy-ion dynamics, Brioni 
1991, Ed. W. Greiner, Acad. Press, New- York, 1992 
G. Welke et al, Phys. Rev. C40(1989)317 
W. Bauer, Phys. Rev. Lett. 61(1988)2534 
E. Suraud et al, Nucl. Phys. A580(1994)323 
E. Suraud et al, Nucl. Phys. A519(1990)171c 



Acknowledgements: The authors thank the french 
german exchange program PROCOPE (grant number 
95073) and Institut Universitaire de France for financial 
support. F. Calvayrac is thanked for providing a version 



4 



